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ABSTRACT 

A compressor-face boundary condition that models the un- 
steady interactions of acoustic and convective velocity dis- 
turbances with a compressor has been implemented into a 
three-dimensional computational fluid dynamics code. Locally 
one-dimensional characteristics along with a small-disturbance 
model are used to compute the acoustic response as a function 
of the local stagger angle and the strength and direction of the 
disturbance. Simulations of the inviscid flow in a straight duct, 
a duct coupled to a compressor, and a supersonic inlet demon- 
strate the behavior of the boundary condition in relation to ex- 
isting boundary conditions. Comparisons with experimental data 
show a large improvement in accuracy over existing boundary 
conditions in the ability to predict the reflected disturbance from 
the interaction of an acoustic disturbance with a compressor. 


NOMENCLATURE 

Roman Letters 

C Chord of the compressor fan blade 

M Mach number 

S Circumferential spacing of compressor blades 

U Desired axial velocity at compressor face 

a Acoustic speed 

e Specific internal energy 

p Static pressure 

r Radial coordinate 

t Time 

w, v, w Cartesian velocity components 


M, W r , V0 

Velocity components in fixed frame 

u, Wr » V0 

Velocity components in blade frame 


Cartesian coordinates 

Greek Letters 

a 

Convective velocity response coefficient 

P 

Acoustic response coefficient 

6 

Small discrete difference 

y 

Ratio of specific heats (y= 1 .4 for air) 

r 

Stagger angle 

p 

Density 

<7 

Local Courant number 

e 

Circumferential cylindrical coordinate 

Subscripts / Superscripts 

a 

Solution point adjacent to boundary 

acoustic 

Component due to acoustic disturbance 

b 

Solution point at boundary 

case 

Location on case 

convective 

Component due to convective disturbance 

cf 

Compressor face 

hub 

Location on hub 

n 

Time level for known state 

n+1 

Time level for new state 

passage 

Conditions within the blade passage 

r 

Radial 

X 

Axial 

9 

Circumferential 
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INTRODUCTION 

The aerodynamic analysis of the flow in an aircraft inlet for a 
turbine engine system may involve the analysis of unsteady flow. 
For example, atmospheric disturbances in pressure or flow angle 
ingested into a supersonic inlet may interact with the compressor 
and reflect acoustic waves which may cause instabilities in the 
shock system and cause unstart (Mayer and Paynter, 1994). 

Acoustic disturbances are changes in pressure which travel 
at the speed of sound relative to the local fluid velocity. They can 
travel upstream, as well as, downstream. Velocity disturbances 
are convected at the local fluid velocity, and so, can only travel 
downstream. Here we are interested in disturbances that travel 
downstream through the duct, interact with the compressor, and 
generate an upstream-travelling acoustic disturbance. 

The analysis of flows in inlets using computational fluid dy- 
namics generally use flow domains in which the compressor face 
is a boundary and some appropriate outflow boundary condition 
is applied. The alternative is to include the compressor or the en- 
tire engine into the flow domain (Suresh et al., 1999). The latter 
approach involves considerably more computational resources 
and is not currently considered practical for design studies. 

Several outflow boundary conditions have been used. Since 
the outflow is subsonic, only one physical quantity needs to 
be imposed. Most boundary conditions have held a quantity 
fixed, such as pressure, corrected mass flow (Mayer and Payn- 
ter, 1994), Mach number (Mayer and Paynter, 1994) (Chung 
and Cole, 1995), or velocity (Mayer and Paynter, 1994)(Payn- 
ter, 1998). However, for unsteady flows, such quantities may not 
be fixed. Experiments have been performed to investigate the un- 
steady nature of disturbances interacting with compressors (Fre- 
und and Sajben, 1997). This has led to an analytic study (Sajben, 
1999). Paynter has investigated a compressor-face boundary con- 
dition that models the unsteady interactions of acoustic (Payn- 
ter, 1997) (Paynter, Clark and Cole, 1998) and convective veloc- 
ity (Paynter, 1998) disturbances with a compressor using locally 
one-dimensional characteristics along with a small-disturbance 
model. The reflected acoustic response is computed as a func- 
tion of the local stagger angle and the strength and direction of 
the disturbance. Paynter et al. (1998) implemented this bound- 
ary condition into two one-dimensional codes and demonstrated 
that it improved the simulation of the response to acoustic distur- 
bances. The current paper addresses the implementation of the 
Paynter small-disturbance compressor-face boundary condition 
into the three-dimensional NPARC code (Chung et al., 1999). 

The Paynter small-disturbance Compressor-face boundary 
condition is summarized and an outline is presented of how 
the boundary condition was implemented. The behavior of the 
boundary condition is examined in comparison with existing out- 
flow boundary conditions and limited experimental data. 


FRAMES OF REFERENCE 

The computational fluid dynamic methods operate in a 
Cartesian (x,y,c) frame-of-reference. For axial flow compres- 
sors, a cylindrical (jc,0,r) frame is useful. At the compressor * 
face, the inlet is assumed to be a constant-area annulus with the 
axial coordinate denoted as jc. The circumferential coordinate 0 is 
coincident with the direction of compressor rotation. The radial • 
coordinate r is perpendicular to the x-6 surface. A blade- fixed 
frame has the cylindrical frame rotating with the fan blade. 

SUMMARY OF THE CASCADE COMPUTATIONS AND 
THE FORMULATION OF A COMPRESSOR MODEL 

Paynter (Paynter, 1997) (Paynter, Clark and Cole, 1998) per- 
formed two-dimensional, time-accurate Euler flow 7 simulations 
with step acoustic and convective velocity disturbances interact- 
ing with a two-dimensional cascade. The responses were exam- 
ined with respect to the stagger angle, solidity, axial Mach num- 
ber, blade loading, blade shape, and disturbance strength. Some 
of the assumptions and major results are summarized here with 
respect to the formulation of a model for the compressor face. 

1 . The compressor face consists of a single row of blades. This 
is equivalent to saying that all reflected waves are due to in- 
teractions with the first stage. Actually, compressors may 
consists of several stages and waves transmitted through the 
first stage may reflect from later stages and propagate up- 
stream through the initial stages. 

2. The volume of the blade passage is negligible relative to the 
volume of the annular section of the inlet. 

3. The rotational speed of the compressor is constant. 

4. The stagger angle T of the compressor blade varies linearly 
in the radial direction, 

r(r) = T hub + r ~ rhuh ( T case -r hub ). ( 1 ) 

r case r hub 

5. u»(v, v e , ve, w', w r , hv ). 

6. The flow properties prior to the disturbance are known. 

7. Acoustic disturbances propagate axially. 

8. Acoustic disturbances yield an acoustic response. 

9. Convective axial velocity disturbances are not possible. 

10. Convective tangential velocity disturbances are possible and 
yield an acoustic response. The effect is a local change in the 
flow angle-of^attack with respect to the compressor blade. 

The response pressure decreases or increases depending on 
whether the velocity relative to the blade is increased or de- 
creased, respectively. 

1 1 . Convected disturbances are convected axially from the “ad- 
jacent” to the “boundary” grid planes without change. * 

12. Radial variations in the velocity do not significantly influ- 
ence the production of the reflected acoustic responses. 
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1 3. The response is not a strong function of blade loading, there- 
fore, the compressor blades are assumed to be unloaded. 
This implies that the flow is attached (unstalled). Further, 
the compressor blades are assumed to be flat and of zero 
thickness. Thus, the flow angle is initially aligned with the 
blades and is equal to the stagger angle. In the blade-fixed 
frame, the tangential velocity component is then given by 

ve C f = u cf tanT (2) 

This allows the passage Mach number to be determined as 
M passage = M c f j COSY (3) 

14. The response becomes insensitive to solidity for solidities 
greater than one. Therefore, the compressor blades are as- 
sumed to have a solidity greater than 1 .0. 

15. If the passage flow is subsonic, the turning of the flow starts 
at the reflected disturbance and ends well into the blade pas- 
sage. If the flow is supersonic all of the turning occurs at the 
reflected disturbance upstream of the blade passage. 

1 6. Four distinct flow regions can be identified and are shown in 
Figs. 1 and 2. Region 1 is prior to any disturbance. Region 2 
is just after the disturbance has passed. Region 3 is between 
the cascade leading edge and the response. Region 4 is in 
the blade passage and downstream of the blade passage that 
is upstream of the transmitted disturbance. 

17. A step increase pressure disturbance causes an increase in 
the axial velocity between regions 1 and 2. The tangential 
component of velocity is unchanged between regions 1 and 
2. For an initially unloaded cascade, this means that the flow 
in region 2 is no longer aligned with the blade passage. For 
weak disturbances, the changes in static flow properties be- 
tween regions 1 and 2 are isentropic. The total properties in 
a blade frame of reference are different in regions 1 and 2 
because the disturbance separates these regions. 

18. For a step change in the tangential velocity being convected 
by the initial axial velocity the axial velocity and static pres- 
sure are unchanged between regions 1 and 2. For an initially 
unloaded cascade, this means that the flow in region 2 is no 
longer aligned with the blade passage. Although the static 
properties are the same in regions I and 2, the total proper- 
ties are different because of the change in velocity. 

19. The tangential component of velocity is unchanged across 
the left running acoustic disturbance. Since only the axial 
component of velocity is changed, the reflected disturbance 
is a plane wave. The changes in static properties across the 
reflected disturbance are isentropic. 



Figure 1. FLOW SCHEMATIC PRIOR TO THE DISTURBANCE / 
BLADE PASSAGE INTERACTION. THE DISTURBANCE APPROACHES 
THE BLADE PAIR FROM THE LEFT. 



Figure 2. FLOW SCHEMATIC AFTER THE DISTURBANCE / BLADE 
PASSAGE INTERACTION. AN ACOUSTIC DISTURBANCE TRAVELS 
UPSTREAM TO THE LEFT. 

ACOUSTIC RESPONSE COEFFICIENT 

The acoustic response coefficient is defined as 

P = (/?3 - Pi ) / (P2 - Pi ) - (4) 

Paynter (Paynter, 1998) and Sajben (Sajben, 1999) have 
each developed an expression for the acoustic response coeffi- 
cient for the response from a compressor face for a small step 
acoustic disturbance. Both approaches derive the response coef- 
ficient by forming a system of linear equations based on the con- 
tinuity, momentum, and energy equations and using the small- 
disturbance and isentropic assumptions along with the character- 
istic relations of the one-dimensional flow. 

Paynter’ s acoustic response coefficient depends on whether 
the Mach number in the blade passage is subsonic or supersonic. 
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For subsonic passage flow (M paS sage < 1), 

W«.r, = .»= (0 ■ (5) 

For supersonic passage flow (M passage > 1), 

mc/j) = i-o. (6) 

Sajben’s acoustic response coefficient is 


SMALL DISTURBANCE BOUNDARY CONDITION 

The two central ideas of the Paynter small-disturbance 
boundary condition are that both acoustic and convective velocity 
disturbances result in upstream-propagating acoustic responses, 
and that the acoustic responses are a simple linear combination 
of the individual responses due to acoustic and convective dis- 
turbances. Thus, a simple time difference can be applied to the 
static pressure at the boundary with respect to the linear combi- 
nation of responses, 

Pb = Pb ^P acoustic + ^P convective' (ii) 

The response due to an acoustic disturbance is 


P(A/ c/ ,r) 


M cf tan 2 r 

\l-M cf ) 1,2(1 +M cf ) + M cf tan 2 T 


(7) 


^Pacoustic — G [ 1 + P ( M c f , F) ] ( p a Pb) • (12) 


One can obtain a non-reflective condition by imposing P = 
0.0. When P = — 1.0, one recovers the constant-pressure bound- 
ary condition. When P = 1.0, one recovers a constant- velocity 
boundary condition. 


CONVECTIVE VELOCITY RESPONSE COEFFICIENT 

The convective velocity response coefficient is defined as 


The subscript a denotes the solution point adjacent to the bound- 
ary in the direction normal to the boundary. The a is the local 
Courant number, 

jfl+1 jTl 

CS - ( u + a) h — — . (13) 

Xb x a 

The time step is chosen such that a < 1 .0 at the boundary. 

The response due to a convective velocity disturbance is 


0 s (P3~P2)/P1 

( v’92 — V0[ ) / v’ei 


( 8 ) 


5p 


convective 


l) 


Pl(^a-^b) 

u c f tanF 4- \% b 


Paynter (1998) has developed a velocity response coefficient 
using similar methods as for the acoustic response coefficient. 
Since axial velocity disturbances result in axial acoustic distur- 
bance, we consider here only step circumferential velocity dis- 
turbances, 5ve- Thus the static pressure, static temperature, and 
axial velocities all remain unchanged. 

The expression for the velocity response coefficient 
a{M c j , F) is dependent on whether the Mach number of the flow 
into the blade passage is subsonic or supersonic. For subsonic 
passage flow (M passage < 1), 


The M x b is the local axial Mach number at the boundary, 

M x b = Ub / <*b- (15) 

Analogous to Eq. 1 1 , the axial velocity can be expressed as, 
“f 1 = 4 + bioacoustic bliconvective (16) 

where. 


a(M c/ ,r) 


Y M?, 

1 -M cf 


tan 



tanf 


(9) 


8u a 


ag[i-P(M c/ ,r)] 

YP£[i+Wc f ,r)] 


8p 


acoustic 


(17) 


For supersonic passage flow (M passage > 1), 


and 


a (M cf ,T) = y M cf (10) 


6 u 


convective 


Y Pi 


5 P 


convective 


(18) 
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IMPLEMENTATION OF THE BOUNDARY CONDITION 

The implementation of the boundary condition into 
the NPARC code required that the conservative variables 
(p, pw. pv\ pvv\ pc) be computed at each grid point on the bound- 
ary at each time step. The computational requirements were no 
greater than required for other outflow boundary conditions in 
NPARC. The steps of the implementation include: 

1. Properties (/?,p,w.v,w 3 a,M) are computed at the “bound- 
ary” and “adjacent” grid points for the n - th time level. 

2. The stagger angle is computed using Eq. I . 

3. The local Courant number is computed using Eq. 13. 

4. The passage Mach number is computed using Eq. 3. 

5. The convective velocity response coefficient a(M c f,F) is 
computed using either Eq. 9 or 10. 

6. The Cartesian velocity components are transformed to the 
fixed cylindrical frame of reference. 

7. The 5 p convective is computed using Eq. 14. 

8. The acoustic response coefficient P is computed using Eq. 5 
or 6. 

9. The acoustic is computed using Eq. 12. 

10. The pressure at the new time level is computed using Eq. 11. 

1 1 . The density at the new time level is computed as 



12. The axial velocity component at the new time level is com- 
puted using Eq. 16. 

13. The cylindrical velocity components in the fixed frame for 
the new time level are computed as 

° (■*.->&>• < 20 > 
and 

° (^tr) - <»)• < 21 > 

14. The velocity components in the fixed cylindrical frame are 
transformed back into the fixed Cartesian frame. 

15. The total energy for the new time level is computed as 

(P^r 1 = (p e)i + ^ (p r 1 - pz) 

+ (p u)l K +I -«?)■ (22) 

16. The procedure is repeated for each boundary grid point. 


RESULTS 

The following results demonstrate the behavior of the Payn- 
ter small-disturbance compressor-face boundary condition in re- 
lation to other boundary conditions with comparisons to the lim- 
ited experimental data for interactions with compressors. 

STEP DISTURBANCES IN A STRAIGHT DUCT 

This case examined the unsteady, inviscid flow in a straight, 
annular duct created when downstream-travelling step distur- 
bances interact with the outflow boundary modelled with various 
outflow boundary conditions. The compressor was assumed to 
have fan blades with stagger angles varying from 45 degrees at 
the hub to 60 degrees at the case. The steady-state compressor- 
face Mach number was 0.2. 

The acoustic disturbance was imposed by initializing the 
flow field in the first quarter of the duct with a +10% step pres- 
sure perturbation. The step propagated downstream as a pure 
acoustic wave. Figure 3 shows the response at a sensor located at 
the three-quarters axial location for the various boundary condi- 
tions. At approximately /=0.62 msec, the incident wave reached 
the sensor and the static pressure steped up 10%. Some “ringing” 
or “oscillation” was noticed, which is due to numerical disper- 
sion. The wave passed by, interacted with the compressor face 
and the reflected, upstream-travelling wave reached the sensor at 
approximately t=\A3 msec. The different behavior of the vari- 
ous boundary conditions became apparent at that time and agreed 
quite well with the expected behavior. The constant-pressure 
boundary condition resulted in an expansion wave that canceled 
out the incident wave. The non-reflective boundary condition 
(p - 0) show that the incident wave exited the outflow bound- 
ary. The time histories for the Chung-Cole and Mayer-Paynter 
constant-Mach number boundary conditions show the correct 
sign for the reflected response with a magnitude approximately 
equal to that of the disturbance. Similar results were obtained 
with the constant- velocity boundary condition (P = 1). The plot 
shows the Paynter small-disturbance boundary condition predict- 
ing a response of about 4% of the steady-state pressure. The plot 
labeled “Sajben” is Paynter’s small-disturbance boundary condi- 
tion using Sajben’s acoustic response coefficient (Eq. 7). 

The radial variation of the acoustic response coefficient P 
for a given compressor-face Mach number can be calculated us- 
ing Eqs. 5 or 7 along with Eq. 1 . if one evaluates the response 
coefficient for the minimum, area-weighted average, and maxi- 
mum stagger angles, one obtains the vertical lines presented at 
the right in Fig. 3. The Paynter expression (Eq. 5) predicts 
a stronger response than the Sajben expression (Eq. 7) and its 
overall response is dominated by the higher stagger angles. 

The responses of the boundary conditions to circumferential 
velocity perturbations of +10% and -10% of the axial velocity 
are presented in Fig. 4. The disturbance reached the sensor at 
a time of / = 3.7 msec; however, nothing was observed since no 
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j 

Paynter Small- Disturbance (-10%) 
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j 

- 
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Figure 3. TIME HISTORY OF THE STATIC PRESSURE FOR A +10% Figure 4. TIME HISTORY OF THE STATIC PRESSURE FOR +/-10% 

STEP PRESSURE PERTURBATION IN A STRAIGHT ANNULAR DUCT STEP CIRCUMFERENTIAL VELOCITY PERTURBATIONS. 


acoustic disturbance is created by the velocity disturbance. At 
a time of t = 6.1 msec the reflected upstream-travelling acous- 
tic wave reached the sensor. Only the Paynter small-disturbance 
boundary condition gave a response. The magnitude of the re- 
flected response was quite small at less than 0.4%, which agreed 
well with the response predicted by Eq. 9. 

UC COMPRESSOR FACE EXPERIMENT 

The experiment by Freund and Sajben (1997) provides the 
limited experimental data for the interaction of an acoustic pulse 
with a compressor. The configuration consisted of a straight, an- 
nular duct connected to a General Electric T58-3 engine by a 
short convergent section. The blades of the first stage of the com- 
pressor had a stagger angle of 45 degrees at the hub and 60 de- 
grees at the case. The Mach number at the compressor face was 
Mach 0.2. A planar acoustic pulse was generated by the rapid 
collapse of a flexible boot located on the hub of the duct. Four 
sensors along the duct between the boot and the compressor mea- 
sured the time history of the static pressure. Sensor 1 is located 
closest to the bump and sensor 4 is closest to the compressor. 

The experiment has been simulated by Paynter, Clark, and 
Cole (1998) using the LAPIN unsteady, quasi-one-dimensional 
flow code along with the Paynter small-disturbance boundary 
condition. Comparisons were made to other outflow boundary 
conditions. Suresh et al. (1999) simulated the experiment by 
coupling NPARC with the ADPAC turbomachinery code. 

The present work simulated the experiment using NPARC 
while assuming unsteady, inviscid flow and using dynamic grids 


to simulate the collapsing bump. The Newton iterative implicit 
method was used with a constant physical time step of 3.0E-06 
seconds, which corresponded to a local Courant number of ap- 
proximately 0.9. 

Figures 5 to 7 show the time histories of the static pressures 
at sensors 1, 2, and 3 for the 7.0 msec span beginning when the 
bump starts to collapse. The lines denoted as “Paynter” used Eq. 
5 for the acoustic response coefficient, P, while the lines denoted 
as “Paynter-Sajben” used Eq. 7. The ranges of values for p are 
indicated in Fig. 5. Also plotted are the experimental data and 
LAPIN results. The LAPIN results were computed using an up- 
dated version which modeled more accurately the bump collapse 
and used a P = 0.362 based on T = 52 degrees and M c f = 0.2067. 
Figure 8 shows the time histories at sensor 2 for the Chung-Cole, 
Mayer-Paynter, and the constant-velocity boundary conditions. 
The time history for the constant-pressure boundary condition is 
not presented, but gave an incorrect response of opposite sign 
with a magnitude equal to the incident pulse. The time history 
for the non-reflective is also not presented, but gave a slight re- 
sponse of the correct sign; however, this response was due to the 
converging duct at the compressor face. 

The plots show the acoustic pulse with a magnitude approx- 
imately 4% of the initial steady-state pressure moving down- 
stream. The initial bump height and collapse time were adjusted 
in the simulations, so that the acoustic pulse matched the ob- 
served pulse. The pulse interacted with the compressor face and 
a reflected acoustic response was generated with an amplitude of 
approximately 1.5% of the initial steady-state pressure. In Fig. 
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Figure 5. THE TIME HISTORIES OF THE PRESSURE AT SENSOR 1 
FOR THE PAYNTER SMALL-DISTURBANCE BOUNDARY CONDITION. 


6, the reflected response is seen at approximately t - 5.3 msec. 
If one sequences through the plots from sensor 3 to sensor 1, 
one can see the reflected response travelling upstream (occuring 
later in time). The Paynter small-disturbance boundary condi- 
tion under-predicts the magnitude of the response; however, the 
reflected response from currently existing boundary conditions, 
as shown in Fig. 8, grossly over-predict the magnitude, which is 
approximately equal to the magnitude of the pulse. The compar- 
isons with the experiment suffer because the simulations do not 
account for the initial guide vanes or later stages of the compres- 
sor which influence the size and shape of the reflected response. 

NASA VARIABLE DIAMETER CENTERBODY INLET 

Simulations of the unsteady, inviscid flow in the NASA vari- 
able diameter centerbody inlet (Wasserbauer, Shaw, Neumann, 
1974) demonstrate the sensitivity of the outflow boundary con- 
ditions in predicting terminal shock motion and tolerance to un- 
start. The inlet is an axisymmetric, mixed-compression inlet de- 
signed for Mach 2.5. A normal terminal shock is located aft 
of the throat. The steady-state flowfield was computed using 
NPARC with the Chung-Cole boundary condition imposing a 
compressor-face Mach number of 0.3 1 . 

The unsteady flow was initiated by a -2% step disturbance 
in the freestream pressure. Figure 9 shows the time history of 
the axial position of the terminal shock. The shock moved for- 
ward as the disturbance passed downstream. The disturbance 
interacted with the compressor face and an upstream-travelling, 
acoustic wave was created, which then interacted with the shock. 



Figure 6. THE TIME HISTORIES OF THE PRESSURE AT SENSOR 2 
FOR THE PAYNTER SMALL-DISTURBANCE BOUNDARY CONDITION. 



Figure 7. THE TIME HISTORIES OF THE PRESSURE AT SENSOR 3 
FOR THE PAYNTER SMALL-DISTURBANCE BOUNDARY CONDITION. 


The constant-pressure boundary condition created a compres- 
sion wave, which results in excessive shock motion with even- 
tual unstart. The Paynter small-disturbance and the Chung-Cole 
constant-Mach number boundary conditions created reflected ex- 
pansion waves as the response. The constant-Mach boundary 
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Figure 8. THE TIME HISTORIES OF THE PRESSURE AT SEN- Figure 9. THE TIME HISTORIES OF THE TERMINAL SHOCK POSI- 

SOR 2 FOR THE CHUNG-COLE, MAYER-PAYNTER, AND CONSTANT- TION SUBJECT TO A FREESTREAM PRESSURE DISTURBANCE. 

VELOCITY BOUNDARY CONDITIONS. 


condition created the stronger response, which results in less for- 
ward motion of the shock. Paynter, Clark, and Cole (1998) ob- 
served this in the one-dimensional analyses and showed that the 
behavior became more apparent with stronger disturbances. The 
result was that the constant-Mach boundary condition would in- 
correctly overestimate the unstart tolerance. 

CONCLUDING REMARKS 

The Paynter small-disturbance compressor-face boundary 
condition was implemented into a three-dimensional computa- 
tional fluid dynamics code and demonstrated a marked improve- 
ment over existing boundary conditions in predicting the re- 
sponse to unsteady acoustic and convective velocity disturbances 
interacting with a compressor face. 
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